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For polarization experiments involving photon counting we introduce a quasi- 
t" , deterministic eigenstate transition model of the analyzer process. Distributions ac- 

cumulated one photon at a time, provide a deterministic explanation for the law of 
Malus. We combine this analyzer model with causal polarization coupling to calculate 
j-^. i photon-photon correlations, one photon pair at a time. The calculated correlations 

exceed the Bell limits and show excellent agreement with the measured correlations 
of [ A. Aspect, P. Grangier and G. Rogers, Phys. Rev. Lett. 49 91 (1982)]. We 
discuss why this model exceeds the Bell type limits. 



I. INTRODUCTION 

Macroscopic electromagnetic field theory leads directly to the law of Malus for light polarization 
measurements. For experiments where the distributions are accumulated from in-sequence photon 

O" 1 

counting, does the law of Malus still hold, and if so, what gives rise to it? Measurements using 



Ok 



a thin polarizing film have supported the Malus formula for in sequence photon counting jlj. In 
these experiments where photons are counted individually in detectors, it is difficult to justify the 
Malus distribution from macroscopic field arguments. Past efforts to explain the law of Malus for in- 
sequence photon counting experiments, as well as photon-photon correlations, with hidden variable 
models have been unsuccessful f|-|J. These models involve a strong probability component. 

Could we successfully explain the law of Malus and measured correlations by replacing the prob- 
ability feature by a deterministic decision process, accumulating the distributions one particle, or 
particle pair at a time? The affirmative answer to this question is the focus of this paper. A de- 
terministic trajectory model || has been previously used to explain "ghost diffraction" patterns 
||. Here, we introduce a quasi-deterministic model of the analyser process. This model involves 
an independent stochastic variable in each analyzer, and a deterministic criterion for selecting the 
polarization eigenchannel in the analyzer. In this model, photons are viewed as field wave packets 
and trajectory calculations are not involved. 

We use this model in detailed calculations to explain the law of Malus for in-sequence photon 
counting. We then use this same model to successfully explain some well known photon-photon cor- 
relations measurements of |7[|J . In [jl0| we use this quasi-deterministic analyzer model to sucessfully 
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explain the proton-proton correlations of as well as the four-angle photon-photon correlations 
of Aspect et al. Q . This model makes use of causal polarization coupling between photons in a pair 
and local stochastic variables. 



II. STOKES REPRESENTATION OF THE EIGENVALUE PROBLEM 



For crystal beam splitters, beam separation is based on eiegenstates, and in particular, the eigen- 
values of the refractive index matrix jl2| . Because of the difference in the eigenvalues for the e and 
o eigenstates, it is possible to separate these states. Because of the dependence on the eigenvalues, 
we call this type of splitter an eigenvalue splitter. The quasi-deterministic model described in de- 
tail below is based on eigcnchannel selection. Stern-Gerlach type analyzers for spinors are likewise 
eigenvalue splitters. In the model described here, we introduce a particular stochastic variable in 
the refractive index matrix. Changing this particular variable does not affect the eigenvalues, so 
that the usual beam separation arguments are unaffected. 

Polarization is conveniently represented by Stokes variables and rotations on the Poincare' sphere. 
Stokes representations for both spinors and vectors are similar and well developed [jl3| . To formulate 
this eigenstate transition model, we first rewrite certain relations contained in the generic two- 
component eigenvalue equation in terms of a convenient Stokes variable representation. All analysis 
here involves a two dimensional complex representation for the field. 

Consider a Hermitian matrix with elements hu = a 7 ^22 = d, hi 2 = /iexp(— i<f>), and h 2 i = 
/iexp(z</>), where a, d, h, and <fi are real. With ^> T — [^xi^yli ^> x = A cos(/3) exp(ta x ) and $f y = 
Asm([3) exp(i(a x + 5), we write the matrix eigenvalue equation h^S = A'S in component form as 
follows. 

(a - A) cos(/3) + h exp(>(<5 - <j))) sin(/3) =0 (1) 



h exp(-i(5 - <j)))cos{(5) + (d - A) sin(/3) = (2) 

In terms of the above variables, the two eigenvalues are given by the following formula. 

A± = (a + d ± y/(a - d) 2 + 4h 2 )/2 (3) 

We employ two Stokes sets S and P that represent the field and matrix respectively. In terms 
of the field variables above, the components of S have the familiar form Sq = ty^ty, Si = 
5oCOs(2/3), S2 — £0 sin(2/3) cos(<5), and S3 = Sq sin(2/3) sin(<5). The components of P are given by 
Pq = y/(a - d) 2 + Ah 2 , Pi = P cos(2o;), P 2 = Pa sin(2a) cos(^) and P 3 = P sin(2a) sin(^), where 
tan(2a) = 2h/(a — d). In terms of these components of S and P, we can extract from equations ([l]), 
(^), and (Q) the following relations. 

P S = ±P S (4) 



S 3 P 2 - S 2 P 3 = 



Si Pi = ±(Pi) 2 VPo 



(5) 
(6) 



Equation (g) means sin(<5 — <j>) — and follows directly from the imaginary component of (Q) and 




The vectors S and P represent points on two Poincare' polarization spheres with common centers 
but with different radii. The radius of the S sphere is given by 5*0 = = A 2 . The radius of the 
P sphere is given by Pq = A + — A_. The Stokes variables rotate with twice the rotation angle of 
the field components. If the space coordinates are rotated by 9 say, the Stokes variables for a spinor 
field are rotated by 6 also, whereas the Stokes variables for a vector field are rotated by 26. 

In this model, an incident field makes a transition to an eigenstate described by equations (^), (|^), 
and (H). From (^) we see that the eigenvalue sign choice in (H) is related to the sign of the product 
Si Pi. We will make use of this relation in the deterministic transition criteria discussed below. 



III. STOCHASTIC ANALYZER VARIABLE 



The traditional e and o rays of classical macroscopic electro-optics correspond to two points on 
opposite sides of the Poincare' sphere indicated by Si = ±So where in the diagonal frame we have 
Pi = Pq for one crystal type. The two eigenvalues, on which the separation decision is made, 
determine the P sphere radius via A + — A_ = Po, but not the direction of P. For a fixed relative 
phase 4>, this gives one free variable for the matrix Stokes vector. 

The point of view here is that the classical matrix Stokes vector used for a macroscopic field of 
many photons does not necessarily represent the matrix Stokes vector experienced by individual 
photons. The first assumption of this model (The distributive assumption) is that the P vectors 
experienced by individual incident pulses are distributed in the one free variable, at least for the 
surface transition region. For a given incident pulse, we randomly select this degree of freedom of 
P from a distribution described below. 

There are two questions related to the depth (z) dependence in this model. First, are the residual 
deviations of P\{z) from the macroscopic field value only surface effects, rapidly decreasing with 
depth into the crystal, or, do they extend throughout the crystal? Both cases give the same correla- 
tions, and both cases give rise to the law of Malus. To correctly describe the correlations, or the law 
of Malus, it is only necessary with the transition criterion described below, that the surface value 
Pi(0) is selected from a distribution described below. For a macroscopic field of many photons, one 
would not expect to have a distributive P vector because of averaging. It is therefore doubtful that 
we can obtain information about this distribution question with experiments using a macroscopic 
field of many photons. 

Second, at what depth into the crystal does the incident field change so that equations (Q), (||), 
and (||) are satisfied? This transition-depth problem in general has been the subject of debated 
over many years p^ |. Unfortunately, the analysis here provides no answer to these questions. The 
presence of a local stochastic variable should have some observable affect. Indeed, it does. Within 
this model, this stochastic variable makes it possible to correctly describe features of both the law 
of Malus and photon-photon correlations. With the distribution selected, the average of Pi(0) is the 
classical macroscopic matrix stokes variable. 



IV. STOCHASTIC VARIABLE SAMPLING 



For a frame attached to the crystal analyzer (say aligned with the e and o rays), we can partition 
the matrix Poincare' sphere with hemispheres indicated by the sign of Pi. However, viewed from a 
frame attached to one analyzer, this partition for the other analyzer is rotated. With respect to an 
analyzers attached frame, we choose P by choosing 2a via 2a = arg where arg = arccos(w) — 7r/2 
and u is selected uniformly from the interval [—1, 1]. For our linear polarizer, we have cos(0) = ±1 
in the crystal frame. Viewed from a frame fixed to the first analyzer, the P' hemisphere axis for 
the second analyzer is rotated from the first by q = 9 for spinors and q = 29 for vectors. In this 
fixed reference frame, sampling for the second analyzer is made using 2a' = arg{u') ± q where u' is 
selected uniformly from the interval [—1, 1], but independent of u. The independence of the selection 
of u and u' represents a local stochastic element of this theory, and is the reason why we call this a 
quasi-deterministic model. 

One should ask if the above sampling distribution is the only distribution that can successfully 
described the data. There are many distributions with which one can exceed the Bell limits in this 
model. The above distribution and a Gaussian distribution with approximately the same width are 
the only two distributions that the author has found to date that correctly describe the data for 
both the law of Malus and the photon-photon correlations. 

V. DETERMINISTIC TRANSITION CRITERION 

The second assumption of this model (The deterministic assumption) is that the incident pulse 
makes a transition to one eigen-channel or the other, and that the choice is made with a deterministic 
criteria based on initial conditions of the incident pulse at the analyzer and the randomly selected P. 
From (H) we see that the sign of S\P\ indicates the eigen-channel choice. We indicate this product 
as a function of the depth (z) into the crystal, as follows. 

T(z) = S 1 (z)P 1 (z) (7) 

The deterministic criteria on which the calculations here are based is that the sign of T(z) after the 
transition is the same as, and determined by the sign of T(0). For a given value for <Si(0) and -Pi(O) 
we make the eigen-channel decision by testing the sign of T(0). In these calculations, we only use 
one of the three equations (^), (||), and (||). Modeling transition criteria using all three equation is 
a subject of ongoing study by the author. 

VI. LAW OF MALUS 

In the study here, we use two deterministic analyzers in sequence, with the second rotated at 
an angle 9 relative to the first. We view the second analyzer in the reference frame of the first. 
The calculations are made one photon at a time. We use the same incident number of photons 
at each relative angle. The notation (+, — ) indicates the configuration for photons that exit the 
(+) channel of the first analyzers and the (— ) channel of the second. To illustrate the statistical 



nature of the distributions, we make calculations for both a small (100), and a large (40000) number 
of incident photons. For the small number case, the accumulated photon counts at each angle for 
the (+, +) configuration are indicated by the solid circles in Fig. |l|. The fluctuations represent the 
approximate Poisson statistics of the sampling. The output for the large number case for the (+, +) 
configuration is indicated by the solid circles in Fig. |[ The solid line represents the scaled Malus 
formula Ncos 2 (0) where N is half the incident photons. The solid squares in Fig. ^ represent the 
accumulated distribution for the (— ,+) case. The dashed line represents the formula Nsin 2 (6). 
From these results, it is clear that this quasi-deterministic hidden variable model clearly gives rise 
to the law of Malus for in-sequence photon counting. 



VII. PHOTON-PHOTON CORRELATIONS 



To make calculations for photon-photon coincidence measurements, we use two deterministic ana- 
lyzer models to represent the two receiving analyzers. Pair counts are accumulated via deterministic 
calculations, one photon pair at a time. The parameters u and u' for the two analyzers are sampled 
as indicated above, but independently of each other. Pair counts for the four different coincidence 

combinations are represented here by N ++ ,iV_| , N |_ and N . For instance is the count 

for the number of pairs with a (+) for the first analyzer an (— ) for the second analyzer. The 

correlations are calculated using the function j(0) = {N ++ + iV__ — — N y)/N fTE] , where 

N = N ++ + AT__ + + N |_ is the total number of pairs counted. For the causal coupling, we 

used Si(0) = Si(0) where the parameter 2f3 was chosen randomly from the interval [0, 2tt\. In the 
correlation calculations here, the deciding information is in the common sign of the first factors in 
the transition functions (0) for each analyzer. It is only necessary that Si(0) and S^O) have the 
same sign. In Fig. || the solid circles represents 2N ++ /N. This curve follows the cos 2 (9) formula, 
and agrees well with the coincidence measurements of [[7]]|], as well as the results reported in || 
for the coincidence data obtained using a narrow spectral bandwidth and a thin crystal for the 
down-conversion process. 

The solid circles in Fig. [| represent the calculated correlation function Jl5|,|| . The four pair-counts 
were accumulated, one photon pair at a time. Calculations are made using about 10000 photon pairs 
at each angle. The solid line in Fig. |] represents the formula of quantum mechanics. The two curves 
show excellent agreement with each other and with the measurements of || . We emphasize that all 
contributions to "f(9) were calculated using the four pair counts only so there is no chance that an 
accidental scale factor can be included. 



VIII. EXCEEDING BELL'S LIMITS 

The causal polarization link between the two photons in a pair produces an angular correlation link 
even with the presence of the local stochastic variables in the analyzers. This angular dependent link 
is why the usual inequality Bell's theorem p^Jl7[ | doesn't apply to this quasi-deterministic analyzer 
model. 



To explain this angular dependence in more detail, consider the transition function T(0) = 
Si(0)Pi(0) = S P cos(2/3) cos(arg{u)) for the first analyzer and T'(0) = Si(0)P{(Q) = 
SqPq cos(2/3) cos(arg(u') + 29) for the second. Both expressions have the same first factor which 
arises because of the causal coupling between the photons. If we don't use this causal coupling, (e.g. 
choose the angles in the first two factors completely independent of each other), the correlations 
vanish. The independent parameters u and v! represents local stochastic variables of this theory. 
Within the distribution range, and depending on the angle 29, there are certain values of v! for 
which the second factor in T"(0) can be negative, and certain values for which this factor is positive. 
Because of these stochastic variables, we cannot know which pair count combination (such as (H — )) 
we have until after the coincident count is made. However, we still have distribution information. 
The number of times the second factor in T'(0) is positive depends on the angle 29. As a conse- 
quence, we have correlations in the distributions which depend on the relative angle. The presence of 
the local stochastic variables, a necessary part of this theory, does not destroy all causal correlation 
information. Wc emphasize that this is a causal non-local effect, and should not be confused with 
instantaneous action at a distance effects. The correlations do not depend on the distance between 
the analyzers, as long as the causal polarization link between the photons remains intact until the 
photons reach the analyzers. 

In comparison with the results here, it is worth commenting on "Bell's theorem without inequal- 
ities" type arguments that have appeared [§,[l8|]. In studies leading to this work, the author has 
calculated correlations with many different hidden variable probability models, and not found one, 
no matter what sampling distribution used, that agrees with the deterministic model and quantum 
mechanics. This finding seems to be consistent with the studies of hidden variable probability mod- 
els by previous authors |§,|l8) . Replacing the deterministic transition criterion at the analyzers with 
probability amounts to throwing away much of the detailed correlations information. 



IX. CONCLUSIONS 



We have introduced a quasi-deterministic analyzer model, and have demonstrated with detailed 
calculations, that the law of Malus and photon- photon correlations can be explained with a causal 
hidden variable theory, via accumulation, one photon or photon pair at a time. If we omit the causal 
coupling, the correlations are lost. If we don't use a stochastic variable in the analyzer we can't 
explain the data for the law of Malus nor the correlations. If we replace the local deterministic deci- 
sion process in the analyzers with one based on probability, the correlations are reduced. In simple 
terms, if we throw away the information, we can't explain the data. Calculations to successfully 
explain some other quantum measurements have been completed and will be reported on elsewhere 




In this paper, we have only considered polarization correlations. The essential features of this 
model are a causal link variable, local analyzer variables, and a deterministic decision criterion. 
Studies using a model with these features to describe energy-time photon correlations |2l| are 
underway. 
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FIG. 1. Malus distribution for few photons. Calculated using the deterministic 
analyzer model: Solid circles indicate the photon accumulation. Counts are 
accumulated one photon at a time. 
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FIG. 2. Malus distribution for many photons. Calculated using the deterministic 
analyzer model: Solid circles indicate the photon accumulation. Counts are 
accumulated one photon at a time. 
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FIG. 3. Normalized photon-photon coincidence counts 2N ++ /N versus relative 
angle: These calculations used the quasi-deterministic analyzer model. The 
distributions are accumulated one photon pair at a time. 
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FIG. 4. Photon-photon correlations calculated using the quasi-deterministic analyzer 
model: Solid circles represent the calculated correlation obtained via accumulation, 
one photon pair at a time. The line indicates the quantum mechanics prediction. 



